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■ Nine classes of integrable boundary conditions for the eight-state supersymmetric model of 

' strongly correlated fermions are presented. The boundary systems are solved by using the coor- 

t-H , dinate Bethe ansatz method and the Bethe ansatz equations for all nine cases are given. 
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^ \ I. INTRODUCTION 

' i Lattice integrable models with open-boundary conditions are one of the recent developments which deserves careful 
elaborations. As many systems in nature are confined in finite boxes (or intervals for one-dimensional systems), the 
effects of boundaries are very significant. This is particularly so for integrable systems since boundary conditions 
generally spoil the integrability of the bulk models. Therefore, the problem of how to extend a bulk integrable model 
to include integrable boundary conditions becomes very important. 

A systematic method for treating integrable lattice models with boundaries, that is the boundary quantum inverse 
scattering method (QISM), has been developed by Sklyanin Q and generalized in J^,^]. Within this framework, the 
integrable boundary conditions are determined by boundary K-matrices obeying the (graded) reflection equations. 

Integrable correlated fermion systems constitute an important class of lattice integrable models, which have recently 
attracted much attention In we proposed two new integrable models with Lie superalgebra <?^(3|1) and 

i quantum superalgebra J7 g [gZ(3|l)] symmetries, respectively. These are eight-state fermion models with correlated 

single-particle and pair hoppings, uncorrelated triple-particle hopping and two- and three-particle on-site interactions. 
. By eight-state, we mean that at a given lattice site j of the length L there are eight possible states: 
> ■ 

|0>, ct + |0), ct |0), ct_j0), 
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where ct Q (c JjQ ) denotes a fermionic creation (annihilation) operator which creates (annihilates) a fermion of species 

■ ol = +, 0, — at site j; these operators satisfy the anti-commutation relations given by {cj a , Cj t p} = 6ij5 a p. 

In a series of papers, we have constructed a large number of integrable boundary conditions for various models of 
strongly correlated electrons In this paper, we are concerned with the integrable eight-state fermion model 

with Lie superalgebra gl(3\l) symmetry. We present nine classes of boundary conditions for this model, all of which 
are shown to be integrable by the graded boundary QISM recently formulated in || . We solve the boundary systems 
by using the coordinate Bethe ansatz method and derive the Bethe ansatz equations for all nine cases. 

This paper is organized as follows. In Section II the boundary model Hamiltonians are described. In the following 
sections we establish the quantum integrability for all these boundary conditions, and derive the corresponding Bethe 
ansatz equations in terms of the coordinate Bethe ansatz method. The last section is devoted to the conclusion. 

X 

^ . II. BOUNDARY MODEL HAMILTONIANS 

We consider the following Hamiltonian with boundary terms 
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H=J2 H if+i + i^ oundary + ff r b ° undary , (II. 1 ) 
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where H,. onn ^ , H^° undaiy are respectively left and right boundary terms whose explicit forms are given below, and 
•ffj'fji is the Hamiltonian density of the eight-state supersymmetric U model || 

#£5+1 0) = ~ ^''j "''^ : -- + h ' C -) ° XP \ ~\ ^ + n i+ l >^ + | 51 ("j"./3 n 3.7 + ™3 + l,/3™j+l, 7 ) 

~ 2(g+ 1) 51 ( C J> C 3,/3 C J + 1 ^ C J + 1 .« + hX ') eXP { _ f ( n i-7 + ^3+l,7)| 

+ ^2( n 3,a + n j+ha ) - — ^("J>"J^ + n J+l,c« n 3+l,/3) 
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( n 3,+ n j,onj~ +nj + i t+ rij + i t oiij+i-), (II. 2) 



(ff + l)(ff + 2) 

where nj a is the number density operator rij a = Cj a Cj a , rij — rij + + iijq + rij_; and 

T? = -ln-^- > C = Hg+l)-iln.g(.g + 2), £ = -ln-£-. (II.3) 
g + 1 2 3 + 2 



We claim that the boundary Hamiltonain ( II. 1 ) is integrable under the boundary conditions: 

Case (i) : ij ; b 4 ° undai y = - - j 3 ^ frn - A( Tll+ „ 10 + mom- + + ^T^~^Tj ni + ni o n i-^ 

^boundary = _ 2 ff | „ L _ 2 ( nL+ni0 + n L0 n L _ + n L+ n L _) + J * n L+ n L0 n L _ ) ; (II.4) 



Case (ii) : i^ oundar y = ^ ^ ni0 + Hl _ _ _?_. rilo7ll _^ 



^boundary = 2^ + _ _J_ nLQnL _ j . , fj - , 

/*''\ j-rboundarv ^5 rrboundarv ^ /tt r*\ 

Case (m) : H H y = jjjj m_, # rt y = — n L _, (II.6) 



plus six others, numbered from Case (iv) to Case (ix) below, which are built from the above three cases by using the 
fact that boundary conditions on the left and on the right end of an open lattice chain can be chosen independently. 
Throughout, £±(a = /, II, III) are some parameters describing the boundary effects. 



III. QUANTUM INTEGRABILITY 

Quantum integrability of the boundary conditions proposed in the previous section can be established by means of 
the (graded) boundary QISM recently formulated in j| . Indeed, the integrability corresponding to the above Case (i) 
has been shown in pOj . We now establish the integrability for the remaining eight cases. We first search for boundary 
K-matrices which satisfy the graded reflection equations: 

1 2 2 1 

Ri2(ui - u 2 ) K- (ui)i?2i(ui + U2) K- (u 2 ) =K- {u 2 )R\ 2 {ui + u 2 ) K- (ui)i?2i(wi - 112), (III.l) 
R^(- Ul + u 2 ) 4* 1 (u 1 )R 12 (-u 1 -u 2 +A) Kf- (u 2 ) 

=K l f 2 (u 2 )R 2 i(- Ul -u 2 + A) Kf 1 ( Ul )i?f 2 12St2 (-7n + u 2 ), (m.2) 
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where R(u) E End(V ® V), with V being 8-dimensional representation of gl(3\l), is the R-matrix of the eight- 
state supersymmetric U model ||, and i?2i(u) = P12-R12 ( u )-f*i2 with P being the graded permutation operator; the 
supertransposition st p (// = 1,2) is only carried out in the fi-th factor superspace of V ® V, whereas ist^ denotes the 
inverse operation of st^. 

It can be checked that there arc three different diagonal boundary K-matrices Q, KL(u), K h {u), Ki n (u), which 
solve the first reflection equation (HI.l): 



KL(u) = - r 



1 



-^diag (Ai (u),BL («), BL (u),Bi («), C£ («), C£ (u), <?£ («), £>£ (u)) , 



ei(2-e)(2+e) 

^- 0) = 777777^77T dia g B -(u), B H [u), B H [u), B H [u),C H {u),C H {u)) 



e(2-e 

tf£"(u) = ^diag(A£ IJ (u),A£"(u),A£"^ , 



(III.3) 
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(III.4) 

The corresponding K-matrices, iT£(u), j^£ f (M ), iCj j (it), can be obtained from the isomorphism of the two reflection 
equations. Indeed, given a solution K0l(u) of (1II.1), then K+(u) defined by 



Kl s \u) =K°L(-u + 2), a = I, II, III, 



(III.5) 



are solutions of (III. 2). The proof follows from some algebraic computations upon substituting ( III.E) ) into ( III.2| ) and 
making use of the properties of the R-matrix . 

Following Sklyanin's arguments |l]] , one may show that the quantity 71 (u) given by 



71 (u) = T(u)K_{u)T-\-u), T{u) = R, 0L {u) • .-J2bi(«), 
satisfies the same relation as K-(u): 

1 22 1 

.R12 (ui - u 2 ) T- {ui)R 2 \{ui + u 2 ) T- (u 2 ) =T- (u 2 )Ri 2 {ui + u 2 ) T- {u\)R 2 i{ui - u 2 ). 

Thus if one defines the boundary transfer matrix t(u) as 

t(u) = str(K + (u)T-(u)) = str (K + (u)T(u)K_(u)T- 1 (-u)) , 



(III.6) 



(III.7) 



(Illi 



then it can be shown || that [t(u\), t(u 2 )] = 0. Since K±(u) can be taken as K±(u), K±(u) and K^ 1 ^), respectively, 
we have nine possible choices of boundary transfer matrices: 



(u) = str (JTJ(u)T(«)ir^(«)r- 1 (-tt)) , a,b = I, II, III, 



(111.9) 



which reflects the fact that the boundary conditions on the left end and on the right end of the open lattice chain are 
independent. 



1 Non-diagonal K-matrices exist but they are beyond the scope of this paper. 



3 



Now it can be shown that Hamiltonians corresponding to all nine boundary conditions are related to the second 
derivative of the boundary transfer matrix T^ a ' b ^(u) (up to an unimportant additive constant) 



H = 2g H^ b \ 
r («,6)"(0) 



L-l 



A(V + 2W) jr{"' J,3+L ' 2" ' 2(V + 2W) 

+2str U a ' + m *\ +str (k\ (0) (tf*) 2 ) 



str K a + (O)Glo 



(111.10) 



where 



V = str Q Kl (0), W = str \K a + (0)H« 



(III.ll) 



with Pjj+i denoting the graded permutation operator acting on the j-th and j + l-th quantum spaces. More precisely 
Case (i), Case (ii) and Case (hi) correspond to H^' T \ H^ 11 ' 11 ^ and H^ 111 ' 111 ^, respectively. We arrange the remaing 
six cases in the following order so that we have the correspondence: Case (iv) <-> H^ 1 ' 11 ^, Case (v) <-> H^ 11 ' 1 ^, Case 
(vi) <-> H^' IH \ Case (vii) <-> H^ nl ^\ Case (viii) <-> H^ u ' nr > and Case (ix) <-> H<- HI > n \ 

Let us also remark that for general boundary parameters £± the boundary terms listed above break the original 
<?^(3|1) symmetry of the bulk model into U(l) x [/(l) x U(l) symmetry (generated by fermion number operators). 



IV. BETHE ANSATZ SOLUTIONS 



Having established the quantum integrability of the boundary model, we now solve it by using the coordinate space 
Bethe ansatz method. Following O|^,|L|l0l , we assume that the eigenfunction of Hamiltonian (II. 2) takes the form 



l*> = V ait ... taN (xi,---,x N )cl. iai ■■■cl NOCN \0), 

{(xj.aj)} 

N 

(kpqi, ■ ■ ■ , k PQN ) exp(i Y k Pj x j)> 

P j=l 



(IV.l) 



where the summation is taken over all permutations and negations of ki , ■ ■ ■ , fcjv , and Q is the permutation of the N 
particles such that 1 < xq\ < ■ ■ ■ < xqn < L. The symbol ep is a sign factor ±1 and changes its sign under each 
'mutation'. Substituting the wavefunction into the eigenvalue equation H\^>) = E\ty), one gets 



■A... J ctj (' ' ' ) kj : kii ' ' ') — ^ij Q^i'i kj )-A... f cti i a : j (' ' ' : ki % kj 7 ' ' ')' 

Aai,— (~ kj ,•■■) — s L (kj ; pi ai )A aii ... {kj, • • •), 
A... t0li (• • • , kj) = s (kj) PLcti )-<4....,a i (' ' ' j kj), 



where Sij(ki,kj) are the two-particle scattering matrices, 



Sij {ki-, kj j 



9{ki) ~9(k 3 ) + icP t3 
6(h) -0(kj) + ic 



(IV.2) 



(IV.3) 



where Pij denotes the operator interchanging the species variables on and ctj, (a.i,otj = + ,0, — ),the rapidities 0(kj) 
are related to the single-particle quasi-momenta kj by 0{k) = ^tan(^) and the dependence on the system parameter 
g is incorporated in the parameter c = 1/g. s L (kj;pi ai ) and s R (kj;pL ai ) are the boundary scattering matrices, 

1 - PiQ.e 1 ^ 



s (kjipion) = 



s R (kj-,p Lai ) = 



1 -PiQ,e" 
1 -PL ai e 



1 - J3L Qi e lfc J 



= 2ife 3 -(L+l) 



(IV.4) 
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with pi^ and pL ai being given by the following formulae, corresponding to ([1.4 - II.6), respectively, 

2ff 



Case (i) : p 1+ = p w = p x _ = pi = -1 



Pl+ = Plo = Pl- = Pl = -1- 



2 — f- ' 

2.g 



Case (ii) : p 1+ = -1, p 1Q = p x _ = -1 



PL+ = -1, Plo =Pl- = -1 + 777; 



Case (hi) : p x+ = p w = -1, pi_ = -1 



Pi+=PiO = -l, PL— = —1 1 /JX . 



2-a 
2ff 

25 
^_ 



1 : 



(IV.5) 



(IV.6) 



(IV.7) 



As is seen above, the two-particle S-matrix ( IV. 3 ) is nothing but the R-matrix of the <?^(3)-invariant Hciscnberg 
isotropic magnetic chain and thus satisfies the quantum Yang-Baxter equation (QYBE), 



Sij {ki, kj*)Su{ki : ki)Sji {kj^ki) — Sji (kj , ki)Sn(ki, ki)Sij iki-, kj*). 
It can be checked that the boundary scattering matrices s L and s R obey the reflection equations: 
Sji(~kj, —ki)s L (kj ; p± aj )5Vj (— fcj, kj)s L (k^; P\ ai ) 

— S iki]Plai)S ji( kj,k{)s {kj\p\ a ^)Sij{ki^kj^ 
Sji( fcj", ki)s (kj\PLaj)Sij(ki, kj)s (ki)PLai) 

= s (ki\ PLai)Sji{kj , —ki)s {kj]PLai)\Pai)Sji(kj^ki). 
This is seen as follows. One introduces the notation 



(IV. 



s(k;p) 



1 — pe 



— ik 



1 — pe 1 



(IV.9) 



(IV.10) 



Then the boundary scattering matrices s L (kj\p\ ai ), s R (kj;pL ai ) can be written as, corresponding to (IV.5 - [V.7] ), 
respectively, 



Case i : s L (kj-,pi ai ) = s(—kj]p\)I , 

s R {k ] -p Lat )^e lk ^ L+1 h(k 3 -p L )I- 

I 

Case ii : s L (k j ;pi olt ) = s(-kj-,pi + ) I (Z-x] 



(IV. 11) 












s H (k j]PLm ) = e^ L+1 h(k riPL+ ) 



/1 



V 

1 

Case hi : s L (kj\p\ ai ) — s(— fc 3 ;pi + ) C 1 







(IV.12) 







K_ — Xj 

1 



^(fc,;pL Qs ) = e ^ 2(i+1) S (%;pL + ) 1 



^±-^ 



(IV.13) 



where J stands for 3x3 identity matrix and pi+, pl+ are the ones given in (IV.6); £±, k± are parameters defined by 

(IV. 14) 



c± = — — — ) K ± ~ 



2fj 



2fj 
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The boundary scattering matrices for cases (iv) - (ix) can be easily built from (IV. 11), (IV. 12) and (IV. 13). We 
immediately see that ( |[V.ll|) are the trivial solutions of the reflection equations ( IV.9| ), whereas (IV. 12) and ( IV.13| ) 
are the diagonal solutions plp| . 

The diagonalization of Hamiltonian (II. 2) reduces to solving the following matrix eigenvalue equation 

Tjt = t, j = l,---,N, 

where t denotes an eigenvector on the space of the spin variables and Tj takes the form 

Tj = Sj (kj )s L (-kj ; via, )RJ (kj)Rj (k 3 )s R (k 3 ■ p L<7j )S+ (kj ) 

with 



(IV. 15) 



(IV.16) 



Sj~(kj) = Sj,N{kj,kff) ■ ■ ■ Sjj+i(kj, fcj+i), 
Sj (^j) = Sj,j-i{kj, kj-i) ■ ■ ■ Sj t i(kj, fei), 
Rj( k j) = Si^ih, -kj) ■ ■ ■ Sj- h j{kj-i,-kj), 



R j ( k j) - Sj+x,j(kj + i,-kj) ■ ■ ■ S N j(k N , -kj). 
This problem can be solved using the algebraic Bethe ansatz method. The Bethe ansatz equations are 

- A^ 1} + ic/2 9j + + ic/2 



N A (1) -0 



e ik ^ L +^F(kj; Pl+ ,p L+ ) 
ic/2 Xi 1] + 6 j + ic/2 



ic/2 + 6 3 - ic/2 



Mi 

JJ B 3 - - ic/2 6, + 

Mi 

G(A«;C-,C + )n 



ic/2 



— A^ 1} - A (1) 



+ ic X^ + \X> + 



(i) 



(i) 



p=i 

p^a 



IC W + Xp' 



ic 



M 2 (1) 

P =i a < 

a = 1 



A 



(2) 



ic/2 X^ } + X p 2) 



ic/2 



(i) 



-A^ 2) + ic/2 X£ } +X { p 2) + ic/2 



AD 



X {2) - X ( P 1] + ic/2 X (2) + A^ + ic/2 
11 A2) ,(D .._,„, (2) , ,(i) ic/2 



M-2 



x {2) - xy + ic xy + xy> + 



p=l A 7 



ic/2 X'-y' + X), 



K (^ ] ;«-.«+) n ,( 2) 

p=i x 1 



v (2) 



(2) 



k (2) 



A 



(2) 



ic X (2) + X 



(2) 



P#7 



,Af 2 



ic 



(IV. 17) 



(IV.18) 



where 



F(kj;p 1+ ,p L+ ) = s(k j ;p 1+ )s(k j ;p L+ ), (for all cases) 
1 case (i) 

— SLZA^t 2 2_L cage ^j Z j 



G(AW;C-,C+) = < 



(C--A<, 1) + f ) (C-H-A^ + f) 
1 



(C++A< r 1 ' + ^) 



(C+ 

(C- 


-A^ + f) 
+Ai 1) + ¥) 


(C- 


-A^+f) 


1 




1 




(C- 


+A< 1) + f) 


(C— A^+f) 

(e + +A«+f) 



^ (C + -AW + f) 



case (m) 
cose (w) 

case (u) 

case (in) 
case (wn) 

case (urn) 
case (ix) 
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1 
1 

(K_+A< f 2) +ic) ( K+ +Ay ) +ic) 

(K_-A<, 2) +ic) (K + -\^ ] +ic) 

1 
1 



(k + -A* 2, +ic) 



O- 



(K++A< r 2) + lc ) 

T2l 

— Ai '+ic) 

(K_+Ay>+ic) 



(«- 



case 
case 

case 

case 
case 



case 



case 



(0 
(it) 

(Hi) 

(iv) 
(v) 

(vi) 
(vii) 
(viii) 
(ix) 



(IV.19) 



The energy eigenvalue E of the model is given by E = —2 J^j=i cos (modular an unimportant additive constant 
coming from the chemical potential term). 



V. CONCLUSION 



In conclusion, we have studied integrable open-boundary conditions for the eight-state supersymmetric U model. 
The quantum integrability of the system follows from the fact that the Hamiltonian may be embedded into a one- 
parameter family of commuting transfer matrices. Moreover, the Bethe ansatz equations are derived by use of the 
coordinate space Bethe ansatz approach. This provides us with a basis for computing the finite size corrections (see, 
e.g. fTl)]) to the low- lying energies in the system, which in turn allow us to use the boundary conformal field theory 
technique to study the critical properties of the boundary. The details will be treated in a separate publication. 
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